A Discontinuous Galerkin Method for Higher-Order Ordinary Differential Equations Applied to the Wave Equation
نویسندگان
چکیده
(ABSTRACT) We propose a new discontinuous finite element method for higher-order initial value problems where the finite element solution exhibits an optimal convergence rate in the L 2 norm. We further show that the q-degree discontinuous solution of a differential equation of order m and its first m − 1 derivatives are strongly O(∆t 2q+2−m) superconvergent at the end of each step. We also establish that the q-degree discontinuous solution is O(∆t q+2) superconvergent at the roots of (q + 1 − m)-degree Jacobi polynomial on each step. Furthermore, we use these results to construct asymptotically correct a posteriori error estimates. Moreover, we design a new discontinuous Galerkin method to solve the wave equation by using a method of lines approach to separate the space and time where we first apply the classical finite element method using p-degree polynomials in space to obtain a system of second-order ordinary differential equations which is solved by our new discontinuous Galerkin method. We provide an error analysis for this new method to show that, on each space-time cell, the discontinuous Galerkin finite element solution is O(∆x p+2) + O(∆t q+2) superconvergent at the tensor product of the shifted roots of the Lobatto polynomials in space and the Jacobi polynomial P 2,0 q−1 (t) in time. Then, we show that the global L 2 error in space and time is O(∆x p+1) + O(∆t q+1) convergent. Furthermore, we are able to construct asymptotically correct a posteriori error estimates for both spatial and temporal components of errors. We validate our theory by presenting several computational results for one, two and three dimensions. Dedication I would like to dedicate this work to my father Bechir, my mother Hajer, my brother Yosri and my wife Zayneb. iii Acknowledgments I would like to declare that I owe my success and achievement to God, without whom I would not be able to reach this level of education. I am grateful to my advisor Dr. Slimane Adjerid for his scientific and moral support. For me Dr. Slimane was more than my advisor. Since our first days in Blacksburg, he provided us, my dear friend Mahboub and I, all the conditions to guarantee our success. More than that, he considered us as member of his family and with him We didn't feel like strangers away from our homes. Words cannot express my love and respect to Dr.Slimane. I would like to …
منابع مشابه
Temporal convergence analysis of a locally implicit discontinuous Galerkin time domain method for electromagnetic wave propagation in dispersive media
This paper is concerned with the approximation of the time domain Maxwell equations in a dispersive propagation medium by a Discontinuous Galerkin Time Domain (DGTD) method. The Debye model is used to describe the dispersive behaviour of the medium. We adapt the locally implicit time integration method from [1] and derive a convergence analysis to prove that the locally implicit DGTD method for...
متن کاملA High Order Approximation of the Two Dimensional Acoustic Wave Equation with Discontinuous Coefficients
This paper concerns with the modeling and construction of a fifth order method for two dimensional acoustic wave equation in heterogenous media. The method is based on a standard discretization of the problem on smooth regions and a nonstandard method for nonsmooth regions. The construction of the nonstandard method is based on the special treatment of the interface using suitable jump conditio...
متن کاملFinite Difference and Discontinuous Galerkin Methods for Wave Equations
Wang, S. 2017. Finite Difference and Discontinuous Galerkin Methods for Wave Equations. Digital Comprehensive Summaries of Uppsala Dissertations from the Faculty of Science and Technology 1522. 53 pp. Uppsala: Acta Universitatis Upsaliensis. ISBN 978-91-554-9927-3. Wave propagation problems can be modeled by partial differential equations. In this thesis, we study wave propagation in fluids and...
متن کاملCombined DG-CG Time Stepping for Wave Equations
The continuous and discontinuous Galerkin time stepping methodologies are combined to develop approximations of second order time derivatives of arbitrary order. This eliminates the doubling of the number of variables that results when a second order problem is written as a first order system. Stability, convergence, and accuracy, of these schemes is established in the context of the wave equat...
متن کاملA Discontinuous Galerkin Method for Ideal Two-Fluid Plasma Equations
Abstract. A discontinuous Galerkin method for the ideal 5 moment two-fluid plasma system is presented. The method uses a second or third order discontinuous Galerkin spatial discretization and a third order TVD Runge-Kutta time stepping scheme. The method is benchmarked against an analytic solution of a dispersive electron acoustic square pulse as well as the two-fluid electromagnetic shock [1]...
متن کاملHigh Order Accurate Runge Kutta Nodal Discontinuous Galerkin Method for Numerical Solution of Linear Convection Equation
This paper deals with a high-order accurate Runge Kutta Discontinuous Galerkin (RKDG) method for the numerical solution of the wave equation, which is one of the simple case of a linear hyperbolic partial differential equation. Nodal DG method is used for a finite element space discretization in ‘x’ by discontinuous approximations. This method combines mainly two key ideas which are based on th...
متن کامل